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1 .0 SUMMARY 

This document presents closed-form solutions for a two-burn orbit insertion, 

AOA, and ATO maneuver targets, and for time of ignition of a one-burn deorbit. 
Based on the assumption that the orbits involved deviate only to the first order 
about a reference circular orbit, these solutions are nearly fuel-minimum. They 
are expressed in terms of linear terminal velocity constraint (LTVC) Shuttle 
guidance targets and, thus, are readily applicable to both ground and onboard 
software use. In addition to potential application to real-time targeting, 
the equations may be used as a mission design aid in preliminary definition 
of target loads and definition of abort mode boundaries. 
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2.0 INTRODUCTION 

Guidance targets for the post-MECO orbit insertion, AOA, and ATO maneuvers are 
currently determined premission and loaded into the onboard and ground computers 
as constants. 

A simple numerical scheme for solving the impulsive two-burn problem, assuming 
tangency of the second burn to the final orbit, is presented in reference 1. 
Using this scheme, which requires iterating until a necessary condition of 
optimality is satisfied, suitable guidance targets are obtained. Still lacking, 
however , is a scheme to determine optimum ignition time for a single deorbit 
burn to cover late AOA aborts. 

Although the scheme presented in reference 1 satisfies the two-burn require- 
ments, a closed-form solution requiring less computations and no iterations is 
presented in this document. A consequence of the derivation, which is obtained 
by assuming near-circular orbit transfers, implies that the second burn of the 
optimum two burns is a horizontal maneuver, which verifies the assumption used 
in reference 1. Additionally, using a similar near circular orbit transfer as- 
sumption, a closed-form solution for optimum ignition of a single deorbit burn 
is derived. 
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SYMBOLS AND ACRONYMS 

abort-once-around 

abort-to-orbit 

Earth gravitational constant 

reference circular orbit radius 

reference circular orbit velocity 

orbit transfer angle 

cotangent of half-0 

sensitivity matrix of state deviations about a reference circular 
orbit 

linear terminal velocity constraint 

initial radial difference divided by r c 

initial horizontal velocity difference divided by V c 

initial radial velocity difference divided by V c 

initial orbit radius at initial location 

initial orbit horizontal velocity at initial location 

initial orbit radial velocity at initial location 

final orbit radius at initial location 

final orbit horizontal velocity at initial location 

final orbit radial velocity at initial location 

radial difference divided by r c after second impulse 

horizontal velocity difference divided by V c after second impulse 

radial velocity difference divided by V c after second impulse 

radial difference divided by r c before second impulse 

horizontal velocity difference divided by V c before second impulse 

radial velocity difference divided by V c before second impulse 

horizontal component of first velocity impulse divided by V c 
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Aoq 

AB 2 

Aa 2 

AV 

AVt 

av 2 

AB 

Aa 

Vh 

r 

ci 

°2 

r iT 

v hiT 

r iT 

r fT 

v hfT 

f. fT 

5pT 

5Bt 

6aT 


radial component of first velocity impulse divided by V c 

horizontal component of second velocity impulse divided by V t 

radial component of second velocity impulse divided by V c 

total velocity impulse divided by V c 

magnitude of first maneuver impulse 

magnitude of second maneuver impulse 

horizontal component of velocity impulse 

radial component of velocity impulse 

horizontal velocity 

radial velocity 

linear terminal velocity constraint target line intercept 

linear terminal velocity constraint target line slope 

initial orbit radius at target location 

initial orbit horizontal velocity at target location 

initial orbit radial velocity at target location 

final orbit radius at target location 

final orbit horizontal velocity at target location 

final orbit radial velocity at target location 

radial difference at target location divided by r c 

horizontal velocity difference at target location divided by V c 

radial velocity difference at target location divided by V c 
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4.0 NEAR-OPTIMUM TWO- BURN SOLUTION 


To derive the closed-form solution for the near-optimum two-burn targets, the 
differences between the final and initial orbits are considered as first-order 
deviations about an arbitrary reference circular orbit. This requires that the 
initial states of both orbits be located at the first-burn ignition point, which 
is fixed in this application. Because the terminal orbits of nominal and ATO 
trajectories are constrained only by desired apsis radius, the final orbit state 
at the first-burn ignition may be determined directly. However, because the ter- 
minal constraints for an AOA are desired radius, velocity and flightpath angle 
(functions of velocity versus flightpath angle target relationship) at a de- 
sired entry angle, these elements must be first mapped from entry interface to 
the first-burn ignition point. 

The approximate state sensitivity matrix that relates state perturbations of 
radius, horizontal and radial velocity after transfer through a central angle 
to initial state perturbations of the same elements about a reference circular 
orbit is given as follows (see appendix A): 


"w 2 + 3 

4 

2W 

W 2 + 1 

W 2 + 1 

W 2 + 1 

-2 

W 2 - 3 

-2W 

W 2 + 1 

W 2 + 1 

W 2 + 1 

2W 

4W 

W 2 - 1 


|_W 2 + 1 W 2 + 1 W 2 + 1 


where W is the cotangent of half of the transfer angle 
Hence , 


(Sp, <58, 6a)T = m’(6Po, $8o, 6a o) T 


( 2 ) 


By defining (<$Po, 58 o, <5ag)T as the difference in the state between the terminal 
and initial orbits at the first impulse , then 

(6p2, 582) 6 a 2)^ = M(6po, 88 q, 5cxq)T • (3) 
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A 


where 6 p 0 = (r f0 - r±o)/r c 

« 0 O = < v hfO ~ V hi 0 )/V c 
<Soi 0 = (r-fo - r i 0 )/V c 
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After the first irapluse is applied, the state differences at the second impulse 
location as a result of the first impulse is impulse 


( 6 P 2 , 5p£, <Sa£)T = M(0, A3 1 , Aa-]) T 


(4) 


By definition, at the second impulse 


<$P2 = <Sp2 

AB 2 = <502 ~ 88£ 

Aa 2 = Saj? - 6a£ 


(5) 


Combining equations (3), (4), and (5) leads to the following: 
( 0 , Ag 2 , Aa 2 ) r = M(<5p 0 , 6g 0 - Ag 1f Aa 0 - Aa-,)? 


( 6 ) 


By defining tan <f > 2 = Aa 2 /AB 2 , and thus Aa 2 = Ag 2 tan <J> 2 , and solving equation 
( 6 ) simultaneously for Ag 2 , Ag 1f and bn* the follow!^ is obtained. 


A 62 = <5p 0 (W 2 + 1)/2(2 -- W tan <J> 2 ) 

A3 -j = -Ag 2 + 6 gp + 6 po 

Aa-| = Ag 2 tan 4>2 + W 6 pg + 8 «o 


The minimum of the sum of the two impulses must be found 
ally expressed as 


This is function- 


AV = AV 1 + AV 2 


Where AV 1 = (Ag^ + Aa^)^ 


( 8 ) 


(9) 
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AV 2 = (Ap| + Aa|)^ 


(10) 


Differentiating equation (8) with respect to tan <|) 2 and w yields: 


dAy _ 2 

d tan 4>2 ( 2 - W tan <j> 2 ) 



AVy 


W + 2 




dAv Aa 2 + W^Po A8 2 
dW (2 - W tan 4> 2 ) ^v 2 


A0 n 

* — + tan §2 
AVt 


1^1 + 

\AVi Av 2 /J 


»• 5P 0 


Aa-j 

A\T[ 


( 11 ) 


( 12 ) 


Equation (8) is minimum when equations (11) and (12) are simultaneously zero. 
By performing the indicated operation, the following necessary condition of 
optimality is obtained. 



(13) 


Except for when 


Aa 


5a 


AV 1 AV 


, the second maneuver is horizontal (i.e., Aa 2 = 0) . 


Results presented in reference 1 indicate that a tangent (near -horizontal) sec- 
ond maneuver is very nearly optimal. • 

Substituting this result (Aa 2 = 0) into equations (7) and (11) yields the 
following: 


A8 2 s 6po (W ^ + 1)/4 
A8i = 68 0 + 5Po (3 - W 2 )/4 
Aa-j = 5aQ + W Spg 


(14) 


and 


A8 1 W - 2 Aa-j = AV-j sgn (6pg) W 


(15) 
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Squaring both sides of equation (15) and factoring yields the following: 


Aa-j (Aa-| W 2 + 4A3 -j W - 4Aai) = 0 


( 16 ) 


Either Ao^ =0, or (Aa-| W 2 + W - 4Aa-|) = ,0 

If Aa-j = 0, then the two- impulse transfer trajectory is near tangent at both 
impulses. This, of course, occurs for the special case in which a gravity turn 
maneuver is optimal. For this case, 


W = -tog/Spo 


(17) 


( Note : Equation (17) is an approximate form of the exact equation derived in 

appendix B for a gravity turn maneuver to achieve a desired apsis magnitude). 

If Aa-|W 2 + 4A3-)W - 4Aa-| = 0, then the solution for W is 


W = 86a 0 / I (468g - 6pg) - sgn(Spg) 


((4<5(3 0 - <$Pq) 2 


(18) 


that was originally derived in reference 2. 

Hence, the location of the second impulse from the first impulse is given as 


0 = IT 


tan" 1 ((2W/(W 2 - 1)) 


To obtain the near-optimum LTVC guidance targets for the first maneuver, the 
state of the final orbit at the first ignition location is extrapolated through 
the angle 9 using Kepler equations of elliptic motion (e.g., equations ( A— 1 ) , 
(A-2) , and (A-3) of appendix A) . The first-burn target radius is defined as the 
final orbit radius at the second-burn location. The target velocity line, 
defined as 


r = c-| + cgVft 
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is chosen to ensure optimality of the resulting solution by defining c-j (target 
line intercept) as the final orbit radial velocity at the second-impulse loca- 
tion and C2 (target line slope) as zero. 

The second-burn ignition time is obtained by first using the LTVC powered flight 
guidance (ref. 3) to determine the finite burn effects of the first maneuver. A 
resulting transfer angle from the predicted cutoff state to the target is then 
determined. Using an appropriate set of Kepler equations for time of flight 
such as presented in reference 4, the second-burn ignition time can be defined. 
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5.0 NEAR-OPTIMAL IGNITION FOR A SINGLE-BURN DEORBI T 

For the analysis of the minimum- impulse solution to the single-impulse deorbit 
problem, first consider the difference between the terminal orbit and the ini- 
tial orbit at the entry interface location. 

6p T = (rf£ - r’i'jO/r'c (19) 

= ( V h f T - V hiT )/V c (20) 

&*X = (ffx - rix)/V c (21) 


For the deorbit problem, the linear terminal velocity constraint is applied to 
express rx as a linear function of V^t as follows: 


rfx = ci + C2VhfT 


( 22 ) 


Thus, equation (21) becomes a linear function of equation (20) as follows: 


6oix = K + C2^8 t 


(23) 


where K = (c-j + C2VhiX - ^iT)/ v c 


The minimum impulse on the initial orbit that satisfies the terminal orbit 
constraints of 6px, <53x> <$ a x> must be located. They are considered as state 
perturbations about a reference circular orbit. Using the approximate state 
sensitivity matrix that relates final state perturbations of radius, horizontal, 
and radial velocity to initial state perturbations of the same elements at 
a central angle 0 about a reference circular orbit,- the following is obtained: 

(<Sp T , <$8 t , k + C2 (S ^t) T = M(0, A3, Aa)T (24) 


where A3, Aa are the impulsive velocity components. 


Solving equation (24) for <S3x> A3, Aa yields the following equations: 


63 T = 2WK + (W 2 + 3)6p T ] /2 (c 2 W - 2) 


(25) 
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A3 = + 5pj 
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( 26 ) 


Aa = - K - C250T + W6p^ 


(27) 

The impulse is defi ,d as follows: 



1 

Av = (A3 2 + Aa 2 )^ 


(28) 

The derivative of AV with respect to W requires the derivative of (SP? with 
respect to W, which is obtained by differentiating equation (25) as follows: 

d63x Aa 

dW " c 2 W - 2 


(29) 

Then, 



dW = C 2 w a - 2 + & + C2K ' {PI 1 /AV 


(30) 

Therefore, for AV to be minimum, Aa must be zero. 



After substituting this result into equation (24), the 
are obtained. 

following equations 


Spx = -1— .i A3 

W 2 + 1 


(31) 

w 2 ~ 0 

63t = J A3 

W 2 + 1 


(32) 

K + c2<S3t = — — — A3 
W 2 + 1 


(33) 
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yjgpor- *' 


After combining equations (31), (32), and (33), a quadratic equation in 
W is obtained as follows: 


1 c 2 W 2 - 2W + 2K/6p T - ^c 2 = 0 


( 34 ) 


Solving this quadratic equation yields the following: 


W = 2 


[l - c 2 (K/6pT - |c 2 )] 


il 


/c 2 


(35) 


Hence, the approximate location of the minimum impulse from entry interface is 
given by 


6 s it + tan" 1 [2W /(W 2 - 1)] (36) 

and the approximate Av is given by 

Av = V c 1 6p T (W 2 + 1)/4 | (37) 


The deorbit maneuver ignition time is obtained by first determining the range 
angle from the current vehicle state to the maneuver location. It is given as 

0 tig = 0 T - 0 

where 9? is the range angle from the current state to the deorbit maneuver tar- 
get location. Then, using the current state and the resulting angle as inputs 
to the Kepler equations for time of flight, the time of ignition may be deter- 
mined. 
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6.0 CONCLUDING REMARKS 

By expressing the initial conditions of the initial and final orbits at a common 
location and assuming that their state differences are deviations about a 
reference circular orbit, simple near-optimum solutions for the guidance 
targets for ascent and abort maneuvers are obtained. This simplicity, in 
addition to the current availability of analytic Jg compensation in the on- 
board' program, make these solutions attractive for targeting and may aid 
in mission design. The way that the analytic J2 compensation scheme can 
be used to bias these targets so that the vehicle can fly through the actual 
targets will be the subject of a forthcoming memorandum. 
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APPENDIX A 

SENSITIVITY MATRIX FOR STATE PERTURBATIONS 
ABOUT A CIRCULAR ORBIT 

The sensitivity matrix derived herein for extrapolating state deviations about 
a circular orbit has been widely used in orbit perturbation analyses. It is 
presented here as an easy reference for the applications of this document. 

The approach taken is to first express the basic Kepler equations for elliptical 
orbits in terms of radius, horizontal, and radial velocities, and cotangent of 
half of the transfer angle; second, differentiate them with respect to the ini- 
tial conditions’, and finally, apply circular orbit conditions to the derivatives 
to obtain the desired result. 

The basic Kepler equations for elliptical orbits may be expressed as follows: 


v h = [(W 2 - 1) v h0 - 2Wf 0 + 2y/r 0 V ho ] /(W 2 + 1) 

(A-1) 

f = [2W(V h0 - y/r 0 V ho ) + (W 2 - 1) r G ] /(W 2 + 1) 

(A-2) 

r = r o v ho/ v h 

(A-3) 

where r 0 , V h0 , r 0 are the initial radius and horizontal and radial velocities 
respectively, and W is the cotangent of one-half of the transfer angle. ^ 

Let r c be an arbitrary reference circular orbit radius. Then 
is the reference circular orbit velocity. Define as follows: 

V c = (y/r 0 ) 2 

8 = v h /v c 


a = r/V G 

(A-4) 

P = r/r c 


Equations (A-1), (A-2) , and (A-3) then become 


8 = [(W 2 - 1)8 0 - zm 0 + 2/P 0 8 0 ] /(W 2 + 1) 

(A-5) 

a = [2W (8 0 - 1/P 0 8 0 ) + (w 2 - 1) a 0 ]/(W 2 + U 

( A— 6 ) 

P = P 0 8 0 /8 

( A— 7 ) 


A-1 
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In order to determine the sensitivity of |3 , a, p to $ 0 , a 0 , p 0 

derivatives of equations ( A— 5 ) , (A-6) , and (A-7) are required. Then, 
in 8, a, p can be functionally related to deviations in 8 0 , a 0 , 

follows: 



dp 

dp 

dp 

dP 0 

d&o 

o 

£ 

d 8_ 

d(3_ 

dB_ 

dpo 

dB 0 

da 0 

da 

da 

da 

dP 0 

d3 0 

t O 
<3 
I 'O 



The required derivatives are expressed as follows: 


d 8_ 

- 2/ (p 0 ^o) 

dPo 

W 2 + 1 


2L - 

2 i 

W - 1 - 2/ 1 

( 2\ 
\Po$o) 

dP 0 

+ 

OJ 


d8_ 

2W 


da 0 

W 2 + 1 

/ 2 \ ” 


da 

2W/ (p 0 8 0 ) 


dP 0 

W 2 + 1 

/ 2 \ 

da 

2W (1+1/1 


d&o 

W 2 + 1 


da 

W 2 - 1 


d a 0 

W 2 + 1 



A- 2 


the 

deviations 
Po as 


(A-8) 


(A-9) 


( A-10 ) 


( A— 11) 


(A-12) 


(A— 13 ) 


(A-14) 
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dp 

dP 0 


(1 - ~ p 0 /8) 8 0 /8 

dP 0 


(A-15) 



a - p 0 /3) p 0 /3 

dp 0 


(A-16) 



2W P 0 B 0 /3 2 
W 2 + 1 


(A-17) 


If the initial conditions of equations (A-1), ( A— 2 ) , and ( A— 3 ) ; i.e., r 0 , 

Vho> and r OJ represent a circular orbit and r c is arbitrarily defined 
as r 0 , then 

P 0 = P = 1 

1 

3o = Po 2 = P = 1 

a 0 = a = 0 (A-20) 


( A— 1 8 ) 
(A-1 9) 


The required derivatives then become 


d3_ 

dP 0 

d8_ 

dP 0 


dB 

da 0 


cP 


2 

W 2 + 1 

W 2 - 3 
W 2 + 1 

2 W 
W 2 + 1 


2 W 
W 2 + 1 


( A— 2 1 ) 


(A-22) 


(A-23) 


(A-24) 
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wjgt. pnr* '■•cm 1 --" — • • 


d°t 

4 W 

dB 0 

W 2 + 1 

da 

W 2 - 1 
W 2 + 1 

dp 

dP 0 

W 2 + 3 
W 2 + 1 

dP 

4 

d8 0 

W 2 + 1 

dp_ 

2W 


W 2 + 1 


* 

Define as M the matrix in equation (A-8), Then 


W 2 + 3 

4 



2W 

W 2 + 1 

W 2 + 

1 


W 2 + 1 

2 

W 2 - 

3 


2W 

W 2 + 1 

W 2 + 

1 


W 2 + 1 

2W 

4W 


w£_ 

- 1 

W 2 + 1 

W 2 + 

1 

w 2 

+ 1 


(A-25) 


(A-26) 


(A-27) 


(A-28) 


(A-29) 


(A-30) 


A-4 


* s 
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is a nonsingular matrix whose inverse is as follows: 


w 2 + 

3 

4 

2W 

w 2 + 

1 

W 2 + 1 

W 2 + 1 

2 


W 2 - 3 

2W 

W 2 + 

1 

W 2 + 1 

W 2 + 1 

2W 


4W 

W 2 - 1 

W 2 + 

1 

W 2 + 1 

W 2 + 1 


Therefore, since 

* 

(<Sp, 68, 6a) T = M(6p 0 j 68 0 , 6a 0 ) T 


A 


( A— 3 1 ) 


* 

( 6 p 0 , 68 0 > 6 a 0 ) T = M - 1 ( 6 p, 53 , Sa) T 


then 
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APPENDIX B 

TRANSFER ANGLE SOLUTION FOR A MINIMUM IMPULSE 
TO A DESIRED APSIS RADIUS 

A solution, originally derived by Tim Brand, CSDL, for the transfer angle from 
impulse to target for a minimum- impulse trajectory to achieve a desired apsis ra- 
dius is developed in this appendix. Application of optimal control theory to 
tiiis problem shows that the minimum impulse control vector is parallel to the ve- 
locity vector at the impulse. Hence, flightpath angle y is unchanged by the 
impulse. In terms of horizontal and radial velocity components after the im- 
pulse 


Vft = V + cos y + = V + cos y 
r + = V + sin y + = Vj£ tan y 


( B— 1 ) 


where the superscript + indicates postimpulse elements. 

The basic transfer equation from impulse to target relates the postimpulse hori- 
zontal velocity to the transfer angle 0 and the target radius rx as follows: 





(W^ + 1) Sp + 1 + tan Y W 


U 

rVft 


(B-2) 


where W = cot - 
2 


and <$p = 1 - r/rx 


Likewise, the velocity at the target Vx is expressed as 


V T 

- (W 2 + 1) 

Sp + 1 

= _H 

( B— 3 ) 

2 

rx Vx 


Since ry^ 

= rx ^x ; 

V T 

r 

(B-4) 



r>x 


Substituting 

this result 

into equations (B-2) and ( B— 3 ) and solving for 

W 

yields 





W = 

tan y /<Sp 



( B— 5 ) 


B-1 
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from which 


0 = TT + tan 


1 [2W/(W 2 - 1)] 


(B-6) 


